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ABSTRACT 

We calculate finite temperature effects on a correlation function in the two dimensional 
supersymmetric nonlinear 0(3) sigma model. The correlation function violates chiral 
symmetry and at zero temperature it has been shown to be a constant, which gives rise to 
a double-valued condensate. Within the bilinear approximation we find an exact result 
in a one-instanton background at finite temperature. In contrast to the result at zero 
temperature we find that the correlation function decays exponentially at large distances. 



1. Introduction 

During the last few years there has been a lot of progress in understanding the dy- 
namics of supersymmetric (SUSY) theories. Although the main achievements have been 
obtained in iV = 2 supersymmetry, where e.g. the low energy effective action can be 
calculated exactly there also exist some remarkable results in N = 1. Among these 
are the dual descriptions of supersymmetric QCD, originally proposed by Seiberg 0. 

Related to this last issue is the calculation of certain instanton induced correlation 
functions in N = 1 SUSY Yang-Mills theory |3j]. The action for this theory possesses 
a global R-symmetry that is conserved classically, but broken by quantum effects. Even 
though the symmetry is anomalous, there still remains a nonanomalous discrete subgroup. 
Some specific correlation functions that violate chiral symmetry, and therefore vanish to 
all orders in perturbation theory, can be calculated in an instanton background for small 
distances and are found to have a constant value, independent of the actual distance. The 
reason for this nonvanishing value is of course that the instantons change the chirality 
by a definite value. Furthermore, by using supersymmetric arguments it was shown 
that this result actually is exact, for all distances. By the cluster decomposition this 
implies a nonvanishing value of the gluino condensate, (AA), and the nonvanishing value 
spontaneously breaks the discrete symmetry down to Z2- As has recently been noted, this 
condensate gives rise to domain wall solutions 

It is a well known feature that if a symmetry is spontaneously broken at zero temper- 
ature, it is often restored above some, possibly high, temperature. This is believed to be 
the case e.g. for the Higgs mechanism in the Electroweak Standard Model || and for the 
spontaneously broken chiral symmetry in QCD [[?]. In view of this, it is clearly interesting 
to study also the finite temperature effects on the spontaneously broken discrete symme- 
tries in supersymmetry, especially when one considers the solid results that already exist 
at T = 0. 

Since the exact results at zero temperature depend crucially on specific properties 
of the supersymmetric theory, and since these relations in general are not valid at finite 
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temperature, one should not expect the T > results to be constrained in the same precise 
manner as at T = 0. However, the underlying structure of these nontrivial theories may 
still be simple enough in order to obtain interesting results. 

In the SUSY SU(2) Yang-Mills theory at finite temperature, some preliminary results 
were recently obtained |§, but no definite, conclusive arguments could be given. In such 
a case, it may prove worthwhile to study the possible features in a somewhat simpler 
context. For this purpose we will use the two dimensional supersymmetric 0(3) cx-model 
as a toy model. 

The SUSY a-model shares many of the properties of the SUSY Yang-Mills theory. 
Apart from being supersymmetric, they also have asymptotic freedom and instanton solu- 
tions in common ||. Furthermore, in this toy model there also exists an exact calculation 
of a condensate []IU|] , that spontaneously breaks the nonanomalous discrete chiral sym- 
metry Z4 — > Z2 (for a comprehensive review of the SUSY a-model, see e.g. [[H]]). Some 



precaution in the analogy at finite temperature is needed however, since by Peierls' argu- 
ment one would expect a theory in only one spatial dimension to restore a broken discrete 
symmetry at any nonzero temperature. Nevertheless, our hope is that this model still 
could serve as a mathematical laboratory and provide some new insights to the possible 
scenarios these discrete symmetries may undergo at finite temperature, e.g. in the SUSY 
Yang-Mills theory. 

The calculation is also interesting from a more technical point of view, since we find 
that the quadratic fluctuations around the instanton solution cancel between the bosons 
and the fermions, at all temperatures. This is a rather surprising result, since the different 
boundary conditions imposed on the bosonic and fermionic field at finite temperature 
naively seems to destroy such a "supersymmetric" cancellation. Although this behavior 
seems to be a specific property in two dimensions (in contrast to e.g. SUSY SU(2) 
Yang-Mills in four dimensions ||), it supports the reasoning that the SUSY theories 
provide simple and constrained, but still nontrivial, examples of doable models at finite 
temperature. 



The paper is organized as follows. In the next section we set up the model and recall 
some of the main results at zero temperature. In section 3 we generalize the necessary 
ingredients to finite temperature, and section 4 is devoted to the calculation of the finite- 
T correlator, with some of the details given in the appendices. Finally, in section 5 we 
give our conclusions, and comment on the results. 



2. Supersymmetric 0(3) a-model at zero temperature. 

The Euclidean action of the supersymmetric 0(3) a- model in two dimensions is defined 
as H 

S = ^J d 2 xd^e^D a <^ a D^ a , (2.1) 

where a = 1,2,3 is the internal isospin index, g the coupling constant, D a the superco- 
variant derivative and $ a a real superfield, 

$ a (f, 9) = ip a {x) + m a (x) + heF a (x) , (2.2) 

satisfying the constraint J2 a ^a{x,6) ^ a {x,d) = 1. This model generalizes the ordinary, 
non-supersymmetric a-model. 

Instead of three variables and one constraint, it is convenient to use a stereographic 
projection and trade the original fields for a complex-valued unconstrained field, 

e = ±-— . 2.3 

1 + $ 3 

The dynamical component fields will then be a complex scalar and its fermionic coun- 
terpart if). The original fields ip a and \P a transform as vectors under 0(3)-rotations, and 
based on these transformation properties one can find the corresponding action on <f> and 
if}. Using the three Euler angles < a , 7 < 2n, < (3 < 71 and defining A = tan(/?/2), 



we have 1C 



c * ^ ia - X (1 + A^e^) 

1 + A0e ia V (1 + X(f)e ia ) 2 ' K ' J 



In order to regulate the infrared divergencies it is assumed that the Euclidean space 
is restricted to a sphere of radius R, and the original flat metric then becomes the confor- 
mally flat one, g^ u = fl^S^, with Q — (1 + (x 2 /4R 2 ))~ 1 . At the end of the calculations 
the flat space limit R —>■ oo is understood to be taken. With this modification the action 
fl2.1| ) becomes 



where x — 1 + 4 ) *4> an d 7 M = cr^, ^ = 1 , 2 (a being the Pauli matrices). 

As is well known, the cr-model possesses nontrivial field configurations that extremize 
the action, the instanton solutions ||. Introducing complex coordinates z = X\ + 1x2 
instead of the Euclidean coordinates x\ and X2, any instanton solution 0i nst is characterized 
by an integer k, the topological charge, given by 



where d z = \{d x - id y ). 

For k > 0, the minimal value of the action is obtained for holomorphic fields, satisfying 
<f) — 0, and to this solution corresponds 4k + 2 real- valued bosonic zero modes and 4k 
fermionic ones. Expanding around this instanton solution, 



S 



4 / d 2 x Q 2 X - 2 l^d^* <9,0 + l -Q 3/2 (>7A^o /: V - (d^o 7 V)7^) 



(2.5) 




(2.6) 



<P = (Pmst(z) + —f=<Pq , iff = -0cl + -^V. 
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the action ( |2.5| ) becomes in the bilinear approximation, 




(2.7) 



where Sq = 4-nk/ g 2 is the classical action and xo = 1 + 0fnst0inst- 



Specializing to k = 1, the most general instanton solution is 



y 



+ c , 



(2.8) 
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where \y\ = p corresponds to the size of the instanton, z to its position and c, together 
with the phase of y, to a rotation of the instanton solution in the original internal isotopic 
space. These three parameters are all complex, giving six real collective coordinates and 
hence six bosonic zero modes. There are four real-valued fermionic zero modes, 

#"(*) = flW^-J' , (2.9) 
\z - ZqY 

with #W (i = 1, 2) a Grassmann parameter and a the spinor index. Note that the fermionic 
zero modes, as they stand, are not normalized. 
Now consider the correlation function, 

n (n) (fi, . . . , x n ) = (0|T{O(xi) . . . 0(£ n ) }|0) , (2.10) 

where O = x~ 2, ?/>(l + 75)^; T stands for time ordering and 75 = 03. Due to conservation 
of chirality, it is clear that Il( n ) receives no perturbative contributions. However, similarly 
to e.g. QCD, there exists an axial-vector current that is conserved classically but broken 
by quantum effects, by a "diangle" anomaly The change in the axial charge Q5 is 
given by 

AQ 5 = 4k , (2.11) 

implying the chiral selection rule n = 2k for the correlator. Hence, in a given topological 
sector there is only one correlation function that does not vanish trivially, but can get 
contributions from nonperturbative, instanton effects. 

Using the chiral selection rule, the only correlator receiving a one-instanton contribu- 
tion is n^ 2 ), and this is the correlation function we will consider. When the action (|2.7|) 
(for k = 1) is used in the path integral, the integrations over the bosonic zero modes are 
traded for the collective coordinates c, y and zq, together with the appropriate Jacobian. 
Similarly there is also a Jacobian associated with the fermionic zero modes. For corre- 
lation functions like n^, the entire contribution is saturated by the zero modes in this 
approximation. Therefore we can integrate out the non-zero modes to get determinants 
of the bosonic and fermionic non-zero eigenvalues. The measure for the integral over the 
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collective coordinates is then given by 

I = e~ 47r/9 o d 2 x d 2 c d 2 y d6^ dflttt d6^ dtf™ (J)et'D B )- 1 (Det'D F )J , (2.12) 

where Db and Dp are the operators in the quadratic fluctuation (|2.7|) for the bosonic 
and fermionic part respectively, with the prime indicating that only the non-zero modes 
should be taken into account, and J is the combined fermion and boson Jacobian for 
the transformation to the collective coordinates. The subscript on g indicates that this 
is the bare coupling constant which will be renormalized. Without the zero modes, the 
super symmetric pairing of bosonic and fermionic degrees of freedom and the degeneracy 
of non-zero eigenvalues (that still holds in the presence of the instanton normally 
imply that the boson and fermion determinants cancel. However, in the presence of an in- 
frared regularization this is not necessarily the case. The determinants are also ultraviolet 
divergent, and an UV-regularization is therefore also needed. Now, by using the trans- 
formation rules (|2.4j ) it is easily seen that n^ 2 ^ is 0(3)-invariant, and since the correlator 
is saturated by the zero mode solutions, whose product is invariant under 0(3) transfor- 
mations, the collective coordinate integration measure I has to share this invariance as 
well. Actually the 0(3) invariance requires the contribution from the determinants and 
the Jacobian to depend on the collective coordinates as (in the limit R — > oo) 

(Det'D B )- 1 (Det'D F )J oc \y\~ 2 (l + |c| 2 )" 2 , (2.13) 

and since the integration measure has to be dimensionless, one finds by dimensional 
arguments 

/ = KM 2 e~^ d 2 x ^— j ^| dO^ d0t(U d 0<2) d 0t(2) ; (2.14) 

(l + |c| 2 ) \y\ 

with K a constant and M 2 an UV cut-off. Of course, to find the correct numerical factor 
K one has to perform the explicit calculation. 

The calculation of the correlation function then gives 

n (2) (f 1? f 2 ) = NA 2 , (2.15) 
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where N is a numerical factor and A is the scale parameter in analogy with QCD, 

A 2 = M 2 e- 4 " /9 ° 2 . (2.16) 

This result for the correlation function is a priori reliable only at small distances, Ax = 
\x\ — X2I — > 0, where n( 2 ) is dominated by small-size instantons, but supersymmetric 
arguments show that neither does there exist any multiloop corrections \TE\, nor can it 



ever be x-dependent |TT|. So the result ( |2.15| ) is actually exact, for all values of Ax, and 
by the cluster decomposition this implies a double-valued vacuum condensate, 

(O) oc ±A , (2.17) 

leaving a discrete Z 2 invariance for the condensate. 



3. Preliminaries at finite temperature. 

In this section we generalize the ingredients necessary for a calculation of the correla- 
tion function IT 2 ) at finite temperature. 

At finite temperature, the Euclidean space is restricted to the strip M. defined by 
Ke(z) = Xi e 9ft and < Im(z) = x 2 < (3 = l/7\ where T is the temperature. Although 
the temporal component of the Euclidean space is compact at finite temperature, the 
spatial part still needs to be regulated. Moreover, any infra-red cutoff has to reduce to 
the T = case when the temperature vanishes. The metric will be taken to be conformally 
flat, gyp = Vt 2 5^ u , with 

n= h S^l (31) 

h{/3, R 2 ) + (/? 2 /tt 2 ) I sinh(7rz//3) | 2 ' 1 J 

satisfying Q(x2 + (3) = ^(^2), and where the function h(f3,R 2 ) satisfies 



lim h(8, R 1 ) = h(0, R 2 ) = AR 2 , (3.2) 

/3— >oo 

to ensure lim^^oo Q = Q . The Euclidean action S at finite temperature is then given by 
the same expression as ( |2.5|) , with the replacements Qq —> Q and ff^ cPx — > Jq dx2 dxi, 
and where bosonic (fermionic) fields are periodic (antiperiodic) under x 2 — * Xi + (3. 
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At finite temperature there still exists an exact instanton solution, which can be de- 
duced from the charge k — 1 instanton by adding an infinite string of such instantons, 
located at = nT -1 = n/3 with identical sizes and rotations, 

00 / 1 \ yn 
<Pinst = y — +c=— coth[7r(z-z )/73]+c (3.3) 

n=-oo \ Z ~ Z - ln P J P 

This describes a periodic k = 1 instanton solution with the same number of collective 
coordinates as at T = 0. Note that when ~Re(z — zq) ~~ ► ±°o the solution becomes 
limR e ( z _2 )_ > -|- oo inst = c± yn/P, implying that when T > the collective coordinate y is 
entangled with the isospin rotations parametrized by c. In other words, already at this 
level it is clear that the rotational degrees of freedom will be more complicated than in 
the T = case. This behavior should be contrasted with the case of Yang-Mills theory in 
four dimensions, where there is no mixing of the collective coordinates at any temperature 
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From the instanton solution ( |3.3| ) we find the bosonic zero modes by taking deriva- 
tives with respect to the collective coordinates, and the fermionic zero modes, satisfying 
antiperiodic boundary conditions, are found from the zero temperature solutions 



„=_oo [z-zq- m/3) \ (3 J smh [n(z - z )/(3\ 

and 

4 2) » = 0^ £ (-1)» 7 V —— = 9<*S* ( V 4) cosh 2 P;-^/f] . (3.5) 

At zero temperature, the Greens functions Il( n ) that are not trivially vanishing but 
receive instanton contributions, are determined by the chiral selection rule. Since we are 
interested in the instanton corrections at finite temperature, it is necessary to establish 
a corresponding selection rule at T > 0. Considering the correlation function at finite 
temperature and denoting by Oi = 0(xi), 

(1X0! ■■■O n )) = Jv^V^j V<f)* V<f> (Oi • • • O n )e~ s , (3.6) 

we now perform a global chiral rotation of 

tp(x) tj){x) = e ia75 V(x) . (3.7) 



Note that any transformation of the fields has to respect the appropriate boundary 
conditions. For a generic chiral parameter a(x), the antiperiodicity of the fermionic 
field requires a(x) to be periodic, which is trivially satisfied in this particular case, 
a(x) = a = constant. The crucial step is then to notice that the relation for the ax- 
ial anomaly remains the same at any temperature; intuitively this is rather clear, since 
the anomaly relation can be viewed as a short distance effect and should therefore be 



independent of the influence from the medium |L9], [2(], Using the Fujikawa method 
||22||, we can write 



{T{O x -"O n )) = Jv^Vij'V ( pV ( p*(0' 1 ---0' n )e- s ' = 

= J Vip Dtp Vcf) V(p* exp 2ina — Aia J d 2 x k 



{O x ---O n )e- s , 

(3i 



where k is the topological charge density, 

k = f d 2 xk , (3.9) 



and the integration is performed over the strip M.. The first term in the exponential 
comes from the rotation of 0\ • • ■ O n and the second from the change in the measure. 



Comparing (|378|) with ( |3.6|) , we get the integrated chiral Ward identity 

(2n-Ak)(T(O l ---O n )} = . (3.10) 

As in the zero temperature case, we expand the action around the instanton solution 
up to quadratic fluctuations. The classical action associated with the instanton solution is 
not affected by the temperature, and thus S = 4:ii/g 2 at any temperature. We again get a 
Jacobian associated with the integrations over the bosonic collective coordinates and also 
a Jacobian from the expansion of the fermionic field in Grassmann coefficients, belonging 
to the zero modes. Furthermore, the bosonic non-zero eigenvalues give a determinant for 
a differential operator defined on the space of periodic functions, and the fermionic non- 
zero modes give a determinant for a differential operator defined on antiperiodic functions. 
The correlation function is then calculated to this order by using the integration measure 
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derived from the semiclassical expansion and replacing the fields in Oi by the bosonic 
instanton solution and the fermionic zero modes. At T = there is a cancellation between 
the bosonic and fermionic contributions to the determinants, that can be understood as 
a consequence of supersymmetry. At T ^ the different boundary conditions for the 
differential operators belonging to the bosonic and fermionic sector seem to remove the 
a priori reason for such a cancellation, but as we show below and in Appendix A, it still 
takes place. This calculation is along the same lines as the corresponding one at zero 



temperature, given in [10, 16]. Following that calculation we first remove the singularities 
in the quadratic fluctuations by making the following rescalings: 



sinh W(z — Zq)/0\ 

X \n*J 

tf> = ^ [ ^) sinh 2 [tt(z - z )/(3] 



(3.11) 



A fP 2 



where 



X 



X \7r 



X° ( ~2 ) l sinh ~ z o)/P}\ 2 



(3.12) 



In the bilinear approximation, we now rewrite the action (|2.7|) in terms of the rescaled 
variables (|3.11|) , 



s = s + j d z x ir 

= S + f d 2 xQ 2 



-^- 2 x^x 2 ^x)~^ + ^-^ 

oz oz ] 



X &zx 



(3.13) 



As shown in Appendix A, by defining 
/ 1 



G9/27T) sinh( 27 r.//5) I , f\ = f , ^W® /0) 
V(/3 2 /vr 2 )sinh 2 (^//3) ' v( VOsmUU-, » 



(3.14) 



and 



R\ P = J d 2 x ^fxfp , R 'kl = J d2% -^f'kf'i , 



(3.15) 



the total Jacobian can be written as 

J = 



1 \ DetR 
DetR' 



(3.16) 
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where R' denotes the part from the fermion Jacobian and R from the boson Jacobian. 

As for the calculation of the determinants Det Dg and Det Dp, we find, by varying 
with respect to the instanton parameters (see Appendix A for details): 



5(lnDet'D B ) = 5 (In DetR) + J d 2 x 5 (In X ) 
5 (In Det' D F ) = 5 (in DetR') + J d 2 xd (In x) 



-dudu In x ~ Trdudn In Q 



7T 



2tt 



-d^lnx 

71 



(3.17) 
(3.18) 



The second term in the square brackets of ( 3.17]) is calculated in the limit where the 
IR-regularization is removed, R — ► oo; 



6 {h I HiOd^lnn) =-51n((l + |c + y| 2 )(l 
where y = yn/P- Hence, up to some numerical factor, 
J (Det' Db) 1 (Det'Dp) = M 2 ^-Lj ^- 



I ~|2^ 



(3.19) 



+ \c + y\ 2 )(l + \c-y\ 2 ) / 

where we have inserted the ultraviolet cut-off M. 

The integration measure at finite temperature can now be written as, 

d 2 c d 2 y 



(3.20) 



I = A 2 [ d 2 x -j- 



(l + \c + y\ 2 ) (l + \c-y\ 2 ) \y\ 
and this measure is invariant under 0(3)-transformations, as it should. 



(3.21) 



4. The correlation function 



In this section we perform the explicit calculation of the spatial correlation function 

n(2 ) f x / / ffisi,Q)(i + 75Msi,Q) ^(o,o)(i + 75MQ,o) \ \ r4 n 

1 u \V x 2 (^i,o) x 2 (o,o) ;/ ' 1 ■ ] 

In the semiclassical approximation, when we replace the fields x an d ip by their classical 
values, we have 

a + rirb = ?m(o<n - + o®(-] cosh y {z " Zo)/(3] ) 

X = XO = 1 + Cst^inst , (4.2) 
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and inserting also the integration measure (|3.21|) from the previous section, the correlator 
becomes 

d 2 c d 2 y 



'>) ldX °(l + \c + y\ 2 )(l + \c-y\ 2 ) \y\ 2 



g(i)gt(Dg(2)gt(2) [^l^sinh^/^l 2 
|sinh[7r(z — zo)/f3}\ 4: \smh.[ir z / j3}\ 4 



1 



NA 2 ( ~] / d 2 x 07 — , 

M 4 |sinh[^//3]| 2 



d 2 



{l + 0* nst inst (z)}{l + (0)} 

c oPy 



sinh[7r(z - z )//3]| 4 |sinh[7rz //9]r Li 1 + C s t</w (*)}{! + inst0inst (0)} 



(4.3) 



where N is a numerical factor. Although not written out explicitly, the instanton field 
<f)- mst also depends on c, zq and y. In Appendix B we show how to perform the c-integration 
by using the Fadeev-Popov method, after which the correlation function becomes 



n( 2 )(z) = A%A 2 



7T 



,2 d2 y 

U Xq 2 

\y\ 



1 - 


1 ~|2 

\y\ 


1 + 


1 ~i2 



X 



1 2/ 1 4 |sinh[7rz//3]|' 



|sinh[7r(,2 — ^ )//3]| 4 |sinh[7rz //3]| 



1 



0inst0inst (0)}|c=0 



{(3 2 |sinh[7r(z - z )//?]| 2 + pV |cosh[7r(z - ^o)//5]| 2 } 
{(3 2 |sinh[vrz //?]| 2 + pV |cosh[vr^//?]| 2 } 



x 



(4.4) 



where K is some new constant and p = \y\. Now, writing z = xo+iyo, putting Re(z) = x\. 
hn(z) = and thus neglecting any time dependence, and defining 

\2 



2-KXi _ 2ttx _ 2iry {pn/P) 
x = , x = —— , y = — , t - 



1 , £ 

- and u = Xn 

(pvr//?) 2 + 1 ° 2 



we can write ( [4.4] ) as (dropping the tildes on the integration variables) 

Kit A 2 . /■<» ^ /-l 



n( 2 )i 



-) J du J dy J dt(l-t 2 )\t\x 
1 

(cosh [it — (x/2)] + t cos?/o) (cosh[u + (x/2)] + t cost/o) 
13 



(4.5) 



The remaining integrals are in principle straightforward, although rather nontrivial. Re- 
ferring the details of the calculation to Appendix C, we find 



n< 2 >(a;i) 



Kn 2 A 2 



1 + 27r:r 1 Tcoth(7r:r 1 T) (l - 2 simi^^T)) + 



+ 2 sinh 2 (7rxiT) In (4 sinh 2 (7rxiT)) 



(4.6) 



where we have assumed that X\ > for simplicity. This result is exact in the one-instanton 
background and the bilinear approximation. 

The correlation function is thus seen to depend on the dimensionless combination x{T 
in a rather complicated way, although it is always a decreasing function of XiT, as shown 
in Fig. 1. 




Figure 1: The correlator]! in units of KA 2 as a function of ixx\T . 



However, for special limits it simplifies considerably. For vanishing temperature or 
distance, 



n (2)( 



Xi) — > at A oc A 
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(4.7) 



Since the limit x\ <C T _1 corresponds to distances much shorter than the average sep- 
aration between the constituents of the medium, Ax mean ~ T _1 , there should be no 
temperature effects. Unsurprisingly, (|4.7| ) is in agreement with the direct T = calcula- 
tion. On the other hand, for asymptotically large values, x{F — ► oo, the correlator falls 
of exponentially: 

rt 2 \ Xl ) Xl ^° 2Kn 2 A 2 (nxiT) e~ 2 ^ T . (4.8) 

Note that the inverse correlation length is given by twice the lowest Matsubara frequency, 
ttT. 

5. Conclusions. 

We have derived, within the semiclassical expansion, an analytical expression for the 
instanton induced, finite temperature correlation function in the SUSY 2d nonlinear a— 
model. For large values of the product x{F we find that the correlator is exponentially 
decaying, and when x{F — > it reduces to the well known T = calculation. We would 
now like to comment on the validity of these expressions and any implications for the 
condensate, (O). 

First of all, we have neglected all extra instanton-antiinstanton pairs and also the 
higher order effects beyond the semiclassical approximation. Hence, the result for the 
correlator can only be accurate to the extent that the coupling is small enough for these 
other effects to be neglectable. Therefore, the relevant question is what sets the scale of 
the running coupling. 

At zero temperature, the only scale available is the distance, so in that case the 
semiclassical approximation is a priori trustworthy only at small distances. But super- 
symmetry guarantees that the correlation function has to be independent of the distance 
and thus makes the short distance calculation valid at all distances. However, at finite 
temperature the supersymmetric arguments are not applicable. This is rather obvious, 
since the correlator now depends on x\ even in the lowest approximation. Nonetheless, we 
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believe the semiclassical approximation to be accurate in the high temperature regime. 
The reason is that when x{T ^> 1, the correlation function is saturated by an instanton 
size of the order p < T _1 . Since the relevant scale in the correlation function should be 
set by the instanton size, it seems reasonable to expect the semiclassical approximation to 
be valid when both x{T 3> 1 and T ^> A. Similarly, when x\T <C 1 the correlator should 
be well described by the semiclassical approximation in the limit of vanishing distance, 
x x -> 0. 

If the above scenario is correct, which seems plausible, the exponential decay of the 
correlator should be reliable at high temperatures. By using the cluster decomposition 
we then arrive at 



and the discrete symmetry is restored. 

Since the theory is defined in one spatial dimension, general arguments based on 
the free energy would imply that the symmetry is restored at any positive temperature. 
This argument is not in conflict with the above semiclassical calculation, although it is 
not easy to justify why the approximations should be qualitatively correct even at low 
temperatures. 

Finally, we would like to connect this model to the full, 4d SUSY Yang-Mills theory, 
although we must emphasize that this is highly speculative at the present stage. However, 
if the results of this toy model at finite temperature has any generalizations to the super- 
symmetric Yang-Mills theory, it would indicate that there is a high temperature phase in 
the 4d theory, where the discrete symmetry is restored. Such a restoration would clearly 
be important in connection with the formation of domain walls || , and a phase transition 
at the temperature of symmetry restoration can in that case not be excluded. 
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Appendix A. Calculation of the Jacobian and the determinants. 

In this appendix we calculate the contribution from the Jacobian and the boson and 
fermion determinant at finite temperature. We will follow the T = approach |10, pl| 
and generalize it to the finite temperature case. 

Starting with the boson Jacobian J#, let a M denote the collective coordinates in the 
following way, a,\ = z , a 2 = y, a 3 = c. The Jacobian is then given by 



J R = DetM {B) , M, ( 5 



n 2 (d<t>\ 



inst 



y inst 



XO V 9a L 



(A.l) 



where the integration is taken over the IR-regulated Euclidean strip. By using the explicit 
form of the instanton solution (13. 31) we can write, 



y inst 



' (3 2 \ 

— sinh 2 (7r(2; - zq)/P) = y 



ii) 



m 



., . sinh 2 (7r(2; - Zq)/0) 

smh(2irz /f3) + cosh(2vr2 //?) smh(2nz/f3) - 2 sinh(2vr2 //?) smh 2 (Trz/f3) 
sinh 2 (7r(2; — z )/f3) = 



dO!2 V 7T 

2tt 



da 3 V 71 " 2 / 

3 2 r l 

sinh 2 (7T2 //3) - - sinh(27T0 O //9) smh(2nz/(3) + cos\i(2tcz /(3) sinh 2 (7rz//3) 

71 z \_ 2 



(A.2) 



If we now define the transposed vector 



#=[1, (P/2n)swh(2nz/f3) , (P 2 /tt 2 ) smh\n z/f3) ] 



(A.3) 



and the following matrix 
/ 



U, 



y 



-(f3/27r)smh(27TZ /[3) 
V (/3 2 /7T 2 )sinh 2 (7rz //3) 





cosh(27rz //3) 
-(f3/ it) smh(2nz //3) 



\ 

(tt//3) sinh(27T2 //?) 
cosh(27rz //5) j 



(A.4) 
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we can write the set of equations (|A.2|) in a compact form as 



i,,s ' 1 ^ sinh 2 (7r(^ - z )/P) = V^U • (A.5) 



da \ir 2 I v K 111 ^ 



Furthermore, by rescaling xo — 1 + ^Lst'/W as 



X = Xo^ |sinh[vr(2; - z )/P}\ 2 , (A.6) 
we find that Mj' = where 

= Jd 2 x^rj p . (A.7) 

The boson Jacobian is thus given by 

J B = (DetU) (DetW) (DetR) . (A.8) 

Note that the definitions of / M and U pv reduce to the well known T = formulas in the 
limit of vanishing temperature Wfi . 

In the fermionic case we have to evaluate 

J F = (DetM (F) ) _1 , (A.9) 



where 



M {F) 1 n 
tJ • Xo 



(d 2 x-4^, (A. 10) 

J Yn 



and ipi are the fermionic zero mode solutions, 

, = V^ 2 cosh(7rQ - z )//3) = y^ 1 

Vl P 2 smh 2 (n(z - z )/f3) ' V2 P sinh(7r(^ - z )/P) ' 1 ' ' 

Performing the same calculations as in the bosonic case give 



'P 



2 N 



^ I ^ 1 sinh 2 (vr(z - zo)//3) = U^f* , (A.12) 



where 

j( cosh^nz/P) \ / / ycosh^o//?) -(yir/P) sinh(7rz //3) 
Ji \(p/ic)fnhh(icz/P)J ' ^ \-(yP/ic)smh(icz /l3) y C osh(nz / P) 



(A.13) 
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Now, by putting 

R m= f**^fkfl, (A-14) 
we get = U'^R'U' , and so the fermion Jacobian becomes 

J F = (DetU'^y 1 (Deti?') _1 (Detf/)" 1 . (A. 15) 



Putting the two contributions Jb and Jf from ( |A.8| ) and ( |A.15| ) together, we obtain the 



total Jacobian J, J = JbJf- By explicitly taking the determinants of U and U we finally 
get the result 

v Is/I y DetR 

Now consider the calculation of the determinants appearing after the gaussian integra- 
tion over the quadratic fluctuations, taken over the non-zero eigenvalues only. Beginning 
with the boson determinant, we have to evaluate 

(Det'Dur 1 = exp [- lnDet'£> B ] , (A. 17) 

where the operator Db after the rescaling is given by 

D B = -4Q- 2 xd z x~ 2 d- z X , (A.18) 

and x is as defined in ( |A.6| ) . This expression for the determinant is of course formal, since 
it is divergent and needs to be regularized. In order to do this we will use the proper time 
method, and define the regularized part of the determinant as |TB 

f 00 dt 



In Det'D B = lim 

<E^0 



- r — (Tre-* Ds -p) + a^' 1 - a ln(e) 



(A.19) 



where p is the number of zero modes, being six real-valued ones in our case. The coeffi- 
cients «i and «o are independent of the instanton parameters, as will be verified in the 
t — > limit, so by making a variation with respect to the parameters we get, 

poo / v 

S(\nDet'D B ) = J dt [Tr {8D B ) e~ WB ) . (A.20) 
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Defining D B = —4% d^x Q dzX \ which satisfies D B Vt X&zX — ^ X^zX Db 
allows us to write 



Tr (8D B )e- tr>B = 2Tr [SQnjc) (D B e 



-tD B _ D B e~ tDB 



(A.21) 



and hence 



5(\nDet'D B ) = -2Tr [6(\nx) ( 



e -tD B _ e -tD B 



oo 
t=0 



(A.22) 



Now, when t — > oo it is evident that only the zero modes contribute; there are six real 
valued zero modes associated with D B and none with D B . Thus 



5 (lnDet'Ds) = -2 / d 2 xd(\nx)P^(x) + |im2Tr [SQnjc) (e 



- tT) B _ e - tD B 



(A.23) 



where P B is the projection operator on the space of bosonic zero modes, Pq = <^</v 
with M an orthonormal basis. The remaining trace can be evaluated by a heat kernel 
expansion of G B (z,z) = (z\e~ tDB \z) and G B (z,z) = (z\e~ tDB \z) . To this end, consider 
lim^-^z G B (z, zi) and let it be represented as 



G B (z, zi) = G B \z, zx) [a (z, z x ) + 01(2;, z x )t + . . .] 



(A.24) 



where G^ B is the free field solution. Since G B has to satisfy periodic boundary conditions, 



it is obtained by summing the T = solution: 



G (o )(z z) _& ~ / Iz-z^-in^T 
{Z > Zl) ~ 4^ A' P At 



(A.25) 



The calculation of the expansion coefficients a and a x in ( |A.24| ) is now straightforward 
TTfl. The value of a (z, z) is fixed by the free field case, a (z,z) = 1, and since we are 



ultimately interested in the limit Z\ — ► z, the result for a±(z, z) follows. We get 



Gr(z, z) 



00 

lim e 

n=— oo 
12 



-n 2 p 2 n 2 /u 



fp_ J. 

Ant Air 



+ —d^lnx 



lim 



tt z 1 , 
— + -d ll d ll ]nx 



00 

1 + 2^ e" 

n=l 



A.Tlt + 4:71 



(A.26) 
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From this equation we see that the expansion is actually independent of the boundary 
conditions in the limit t — *■ 0, since the limit t — > is equivalent to /3 — > oo. This result 
can be obtained by a calculation in momentum space as well HI§| . A similar calculation 
for Db gives 



Vt 2 1 1 



(A.27) 



and ( A.23 ) is then given by 

5QnDet'D B ) = -2 / rf 2 x5(lnx)P S (^) + / d 2 x<5(lnx) 



In x - 7q-<V^ m fi 

7T Z7T 



(A.28) 



Since the variation of the determinant is unaffected by the boundary conditions, the formal 
dependence on \ and Q is the same as at T = 0, but both these functions themselves are 
of course very different at finite temperature. 

The remaining part in the boson determinant is to calculate the projection on the space 
of the zero modes. For this purpose, we take as a non-orthonormal basis the periodic 
functions = x -1 //^ The projection operator is then given by P S = Y*n,v dtd^Jlu^ 2 '■> 



with R^v defined in ( |A.7| ). Substituting this into the first term in ( A.28 ) 



-2 J d 2 xd(\nx)P B (x) = -2 J d 2 x 8 {In x) gig 
-2 J d 2 x5(x)T 3 n 2 f:UR-l = 5(\nBetR) . 



(A.29) 



The boson determinant thus becomes 



5 (In Det'D B ) = 5 (In DetR) + J d 2 x 5 (In x) 



In x ~ Trdudn In Q 



(A.30) 



_7T ' ' " 27T 

The fermionic part is calculated much in the same way; after taking the determinant 
over the spinor indices the fermion operator with positive definite eigenvalues becomes 
Dp = — 4n _3 / 2 x^x _2 ^ _1 <9zX^ 1 ^ 2 ) an d varying the regularized expression with respect to 
the instanton parameters we are left with 



5 (In Det'-Djr) = -2Tr [5 (lux) (e 
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-tD F _ e -tD F 



t=0 



(A.31) 



where we have defined Dp = — 4Q 3 ^ 2 x l d z x 2 Q l d z x 1 ^ 1//2 , satisfying 



(A.32) 



As in the bosonic case there is only a contribution from the zero modes in the t — > oo 
limit, and when t —>■ we can make a heat kernel expansion. Expanding in powers of t 
around the free field solution, that now satisfies anti-periodic boundary conditions, we 
find 



lim(z|e Wf \z) = lim 



^ + ±( W n* + iln n ) 



l+2£(-l)"e 



n -n 2 (3 2 n 2 /4t 



n=l 



(A.33) 



The result for Dp is obtained by making the substitution x ~~ > X 



lim(x|e * Df |x) = lim 
t^o t->o 



l+2£(-l)"e 



n-n 2 (3 2 Cl 2 /4t 



n=l 



(A.34) 



Thus we see again that in the limit of vanishing t, the expansion becomes independent of 
the boundary conditions. Including the projection over the zero modes as well, 



5(\nDet'D F ) =-2 / d 2 x 5 (In*) P F (x) + / rf 2 x5(lnx) 



1 



7T 



In x 



(A.35) 



Taking the anti-periodic functions g i = (l/vQ)x 1 f[ as a basis, we can write Pq 
J2i,j g'/diRij 1 ^ 2 , with R'^ being the same as in ( |A.14| ). Hence 



-2 J d 2 x5(lnx)P F (x) = -2 J d 2 x6 (In*) g/g'XjW 

-2 J d 2 x 5{x)T* ff fiR'ifV = & (in Dettf) , (A.36) 



and the variation of the fermion determinant finally becomes 



5 (In Bet' D F ) = 5 (in Detfl') + / d 2 x 5 (In £) 



-dpdp In x 

7T 



(A.37) 
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Appendix B. Integrating out the coordinate c. 



Here we derive the formula (fOj), by using the Fadeev-Popov method. The starting 
point is the expression for the correlation function with all integrations over the collective 
coordinates remaining, 

n (2) W = / , ,2/, f C ^f ( ? , Z^ f^c,y,z ) , (B.l) 
J \y\ 1 + \c + y\ 1 + \c-y\ 



where y = yir / (3 and we have defined a function f(z,c,y,z ) that contains the scale 
factor A, the contribution from the fields in the correlation function and some irrelevant 
numerical factor. 

The basic idea is that by an 0(3) transformation we can always make c = 0. The 
integration over c can thus be replaced by an integral over 0(3) parameters and a suitable 
Jacobian. To find the Jacobian we define 

1 



dQ5{c') , (B.2) 



F(c,y) 

where c' is the value of c after an 0(3)-rotation, 

, = e j T / (c + y)e iQ -X (c-y)e^-X \ 

2 \l + X(c + y)e ia 1 + A(c - y)e ia J ' 1 ' ' 

and dfl is the invariant group measure, which in terms of our 0(3)-parameters reads 

//•2ir pit rco d\ \ 

This construction ensures that F is rotationally invariant. Inserting the trivial factor 

F(c,y) J dQ5(c') = 1 (B.5) 

into the correlation function and performing a rotation to new values c — *> c' , y —> y' and 
zq — > z'q we obtain, by making use of the 0(3) invariance of f(z, c, y, z ), 

n«(*) = [ dQ [ ■ - 2 , ■ d * C ' d y^f° , , -^ F(c',y')f(z,c',y',z' )5(c) 



\yf (l + \c' + yf) (l + \c'-y>\ 2 



^ 2 I . a ^ yd2 ^ a>a -/--(0./y)/(^.0.//..:o). (D.G) 



MM 1 + \y\"~ 
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(JdQ = 2vr 2 ). According to the definition ( gj) , 



lim 



F(0,y) ^oF(c,y) 

r r°° i 

= lim dad-y — ,,. 
c^oJ 'Jo (1 + A 2 ) 2 



dXX 



e i7 ( (c + y)e ia - A (c - y)e ia - A 



f dadXX 
nm27r / — — 

c-o J (1 + A 2 ) 2 



1 - AY e 2 "* 



2 \l + X(c + y)e ia l + X(c-y)e ic 
~6 (c - X{e~ ia + fe ioL )) . 



The delta function can be rewritten as 
S(c-X(e- ia + fe l 
2 



■61 X 



\e~ ia + y 2 e ia \ 



~2 e ia _ e -ia y 2 e -ia _ ^ 
y2 e ia _|_ g-ia ^ 2 e -i« _|_ & ia 

and by inserting this back into ( |B.7| ), 

1 



5 (a — ot(arg(c))) 



2tt 



~i4 1 



^(0,2/) 

With the help of ( p.9|) , we finally write ( [B.6| ) as 

d 2 y d 2 z 



Z = 7T 



1 - 


1 ~|2 


1 4 


1 ~|2 



/(z,0,y,zb) 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



Note that for (3 — > oo, the effect of the c-integration reduces to just an overall constant 
factor, as in UTol. 



Appendix C. Calculation of the correlation function. 



In this appendix we will explicitly calculate the integrals that remain in the expression 
3D for the correlator. We start from 



n< 2 >. 



KnA 2 



sinh 2 (-) J du J dy J dt(l-t 2 )\t\x 



x 



1 



(cosh[w — (x/2)] + tcosyo) (cosh[it + (x/2)] + tcosyo) 



(C.l) 
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where x = 2nx\/ f3. Next we combine the denominators by using the Feynman parameter 
trick and also the addition theorem for the hyperbolic functions. The result is 



nW( Xl ) = 3KttA 2 sinh 2 dt (1 — t 2 )t d<\ (I — <\)<\ x 

/OO r2iv "j 
du I dyo j 
-oo jo (T cosh u + t cos yo ) 



KnA 2 



sinh 2 (^j dt (1 - t 2 )t £ da (1 - a)a x 

roo r2n roo 

I du / dy / rf 77 3 e -7(T cosh cos W ) 
-co JO JO 

= 2iT7r 2 A 2 sinh 2 (^j £ dt(l-t 2 )t £ da(l-a)a d^^ 3 K (^T)I (-ft) , 

(C.2) 



where 



T = ^/cosh 2 (|) - (1 - 2af sinh 2 (£) , (C.3) 

and Kq and Jo are the modified Bessel functions of zero order. The integration over 7 
gives a hypergeometric function, 

J o d 1 ^K ( 1 T)I ( 1 t) = —F(2;2- 1 l;t 2 /T 2 ), (C.4) 



and using the explicit form of the hypergeometric function 

1+^ 2 
(1-z 2 ) 



F(2-2-l-z 2 )=^-^ (C.5) 



the result is, 

IlWfa) = SKn 2 A 2 sinh 2 (|) £ dt (1 - t 2 )t £ da (1 - ct)a ** + * 3 . (C.6) 

Performing the last two integrations, and assuming xi > for simplicity, we finally arrive 
at 

n (2) (jj!) = — ^ — [l + 27ra; 1 Tcoth(7ra; 1 T) (l - 2 sinh 2 (7rxiT)) + 

+ 2sinh 2 (7ra; 1 T)ln(4sinh 2 (7rxiT))] (C.7) 
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